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1. Introduction 

The production of heavy quarks, both hadropro- 
duction and leptoproduction, has become an impor- 
tant theoretical and phenomenological issue. While 
the hadroproduction mode is of direct interest to this 
workshop, [ [lj we shall find that the simpler leptopro- 
duction process can provide important insights into the 
fundamental production mechanisms. [||, |[ f|, || There- 
fore, in preparation for the Tevatron Run II, we must 
consider information from a variety of sources includ- 
ing charm and bottom production at fixed-target and 
collider lepton and hadron facilities. 

For example, the charm contribution to the total 
structure function F2 at HERA, is sizeable, up to 
~ 25% in the small x region. [ [5J Therefore a proper 
description of charm-quark production is required for 
a global analysis of structure function data, and hence 
a precise extraction of the parton densities in the pro- 
ton. These elements are important for addressing a 
variety of issues at the Tevatron. 

In addition to the studies investigated at the Run II 
workshop series]^] we want to call attention to the 
extensive work done in the Standard Model Physics 
(and more) at the LHC Workshop organized by 
Guido Altarelli, Daniel Denegri, Daniel Froidevaux, 



* Contribution to the Physics at Run II Workshops: QCD and 
Weak Boson Physics, Fermilab, 1999. 

* Sub-group convenor. 

In particular, in the Run II B-Physics workshop, the stud- 
ies of Working Group 4 : Production, Fragmentation, Spec- 
troscopy, organized by Eric Braaten, Keith Ellis, Eric Lae- 
nen, William Trischuk, Rick Van Kooten, and Scott Menary, 
addressed many issues of direct interest to this subgroup. 
The report is in progress, and the web page is located at: 
http:/ /www-theory.fnal.gov/people/ligeti/Brun2/ 



Michelangelo Mangano, Tatsuya Nakada which was 
held at CERN during the same period.^ In particular, 
the investigations of the LHC b-production group (con- 
venors: Paolo Nason, Giovanni Ridolfi, Olivier Schnei- 
der, Giuseppe Tartarelli, Vikas Pratibha) and the QCD 
group (convenors: Stefano Catani, Davison Soper, W. 
James Stirling, Stefan Tapprogge, Michael Dittmar) 
are directly relevant to the material discussed here. 
Furthermore, our report limits its scope to the issues 
discussed within the Run II workshop; for a recent 
comprehensive review, see Ref. [ ^| . 

2. Schemes for Heavy Quark Production 

Heavy quark production also provides an important 
theoretical challenge as the presence of the heavy quark 
mass, M, introduces a new scale into the problem. The 
heavy quark mass scale, M, in addition to the charac- 
teristic energy scale of the process (which we will label 
here generically as E), will require a different orga- 
nization of the perturbation series depending on the 
relative magnitudes of M and E. We find there are 
essentially two cases to consider]^ 

1. For the case of E ~ M, heavy-quark production 
is calculated in the so-called fixed flavor number 
(FFN) scheme from hard processes initiated by 
light quarks (u, d, s) and gluons, where all effects 
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http: / /home. cern.ch/~mlm/lhc99/lhcworkshop. html 
J We emphasize that the choice of a prescription for dealing with 
quark masses in the hard scattering coefficients for deeply in- 
elastic scattering is a separate issue from the choice of definition 
of the parton distribution functions. For all of the prescriptions 
discussed here, one uses the standard MS definition of parton 
distributions. 
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of the charm quark are contained in the pertur- 
bative coefficient functions. The FFN scheme in- 
corporates the correct threshold behavior, but for 
large scales, E 3> M, the coefficient functions in 
the FFN scheme at higher orders in a s contain 
potentially large logarithms In" (E 2 /M 2 ), which 
may need to be resummed. [ |?], ||, ||, [To) 

2. For the case of £ > M, it is necessary to include 
the heavy quark as an active parton in the pro- 
ton. This serves to resum the potentially large 
logarithms \n n (E 2 /M 2 ) discussed above. The 
simplest approach incorporating this idea is the 
so-called zero mass variable flavor number (ZM- 
VFN) scheme, where heavy quarks are treated as 
infinitely massive below some scale E ~ M and 
massless above this threshold. This prescription 
has been used in global fits for many years, but 
it has an error of 0(M 2 / E 2 ) and is not suited 
for quantitative analyses unless E 3> M. 

While the extreme limits E > M and E ~ M 
are straightforward, much of the experimental data 
lie in the intermediate region As such, the correct 
PQCD formulation of heavy quark production, capa- 
ble of spanning the full energy range, must incorporate 
the physics of both the FFN scheme and the ZM-VFN 
scheme. Considerable effort has been made to devise a 
prescription for heavy-flavor production that interpo- 
lates between the FFN scheme close to threshold and 
the ZM-VFN scheme at large E. 

The generalized VFN scheme includes the heavy 
quark as an active parton flavor and involves matching 
between the FFN scheme with three active flavors and 
a four-flavor prescription with non-zero heavy-quark 
mass. It employs the fact that the mass singulari- 
ties associated with the heavy-quark mass can be re- 
summed into the parton distributions without taking 
the limit M — + in the short-distance coefficient func- 
tions, as done in the ZM-VFN scheme. This is precisely 
the underlying idea of the Aivazis-Collins-Olness- 
Tung (ACOT) ACOT scheme [ [El] which is based on 
the renormalization method of Collins-Wilczek-Zee 
(CWZ).[ |l^| The order-by-order procedure to imple- 
ment this approach has now been systematically es- 
tablished to all orders in PQCD by Collins. [ 

Recently, additional implementations of VFN 
schemes have been defined in the literature. While 
these schemes all agree in principle on the result 
summed to all orders of perturbation theory, the way 
of ordering the perturbative expansion is not unique 
and the results differ at finite order in perturbation 
theory. The Thorne-Roberts (TR) [ Q prescription 
has been used in the MRST recent global analyses of 
parton distributions. [ The BMSN and CNS pre- 



scriptions have made use of the 0(a 2 ) calculations by 
Smith, van Neerven, and collaborators [ ||, ^) to carry 
these ideas to higher order. The boundary conditions 
on the PDF's at the flavor threshold become more com- 
plicated at this order; in particular, the PDF's are no 
longer continuous across the N to N+l flavor threshold. 
Buza et al.,[^ have computed the matching conditions, 
and this has been implemented in an evolution pro- 
gram by CSN.[ § More recently, a Simplified- ACOT 
(SACOT) scheme inspired by the prescription advo- 
cated by Collins [ 13 was introduced; [ |l6) we describe 



this new scheme in Sec. 
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Figure 1. F| for x = 0.01 as a function of Q 2 in GeV 
for two choices of n as obtained within the 0(a\) FFN 
and (ACOT) VFN schemes. For details, see Ref. [0. 



To compare the features of the FFN scheme with 
the ACOT VFN scheme^ concretely, we will take the 
example of heavy quark production in DIS; the fea- 
tures we extract from this example are directly ap- 
plicable to the hadroproduction case relevant for the 
Tevatron Run II. One measure we have of estimating 
the uncertainty of a calculated quantity is to examine 
the variation of the renormalization and factorization 
scale dependence. While this method can only provide 
a lower bound on the uncertainty, it is a useful tool. 

In Fig. [j], we display the component of F% for the 
s + W — > c sub-process at x = 0.01 plotted vs. Q 2 . 
We gauge the scale uncertainty by varying fi from 
1/2/xoto 2.0 fio with /xq = \J Q 2 + m 2 . In this figure, 
both schemes are applied to 0(a\). We observe that 
the FFN scheme is narrower at low Q, and increases 
slightly at larger Q. This behavior is reasonable given 



4 In this section we shall use the ACOT VFN scheme for this il- 
lustration. The conclusions extracted in comparison to the FFN 
scheme are largely independent of which VFN scheme are used. 
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that we expect this scheme to work best in the thresh- 
old region, but to decrease in accuracy as the unrc- 
summed logs of In" (Q 2 / 'm 2 c ) increase. 

Conversely, the ACOT VFN scheme has quite the 
opposite behavior. At low Q, this calculation displays 
mild scale uncertainty, but at large Q this uncertainty 
is significantly reduced. This is an indication that the 
resummation of the \n n (Q 2 /m 2 ) terms via the heavy 
quark PDF serves to decrease the scale uncertainty at a 
given order of perturbation theory. While these general 
results were to be expected, what is surprising is the 
magnitude of the scale variation. Even in the threshold 
region where Q ~ m c we find that the VFN scheme is 
comparable or better than the FFN scheme. 

At present, the FFN scheme has been calculated to 
one further order in perturbation theory, 0(a 2 ). While 
the higher order terms do serve to reduce the scale 
uncertainty, it is only at the lowest values of Q that the 
0{a 2 ) FFN band is smaller than the 0{a\) VFN band. 
Recently, 0(a 2 ) calculations in the VFN scheme have 
been performed; [ |i| it would be interesting to extend 
such comparisons to these new calculations. 

Let us also take this opportunity to clarify a mis- 
conception that has occasionally appeared in the liter- 
ature. The VFN scheme is not required to reduce to 
the FFN scheme at Q = m c . While it is true that the 
VFN scheme does have the FFN scheme as a limit, this 
matching depends on the definitions of the PDF's, and 
the choice of the fi scale.^j In this particular example, 
even at Q = m c , the resummed logs in the heavy quark 
PDF can yield a non-zero contribution which help to 
stabilize the scale dependence of the VFN scheme re- 
sult^] 

The upshot is that even in the threshold region, the 
resummation of the logarithms via the heavy quark 
PDF's can help the stability of the theory. 

4. Simplified ACOT (SACOT) prescription 

We investigate a modification of the ACOT scheme 
inspired by the prescription advocated by Collins. [ |l3|| 
This prescription has the advantage of being easy to 
state, and allowing relatively simple calculations. Such 
simplicity could be crucial for going beyond one loop 
order in calculations.^ 

Simplified ACOT (SACOT) prescription. 
Set Mh to zero in the calculation of the 
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SACOT- ■ ■ ■ 




5 The gen eral renormalization scheme is laid out in the CWZ 
paper[ hjL The matc hing of the PDF's at 0{a\) was computed 
in Ref. and Ref. [ M. The O(a^) boundary conditions were 
computed i n R ef. [ ij. 
6 Cf., Ref. [u7 for a detailed discussion. 

' bec Ref. r@ for a detailed definition, discussion, and 
comparisons. 



Figure 2. F% as a function of Q 2 in GeV computed 
to 0{a\) in the ZM-VFN, FFN, ACOT, and SACOT 
schemes using CTEQ4M PDF's. Fig. a) x = 0.1 , and 
Fig. b) x = 0.001. Figures taken from Ref. [ @. 

hard scattering partonic functions a for in- 
coming heavy quarks. 

For example, this scheme tremendously simplifies 
the calculation of the neutral current structure func- 
tion F^harm even a ^ 0(a\). In other prescriptions, the 
tree process 7 + c — > c + g and the one loop process 
7+c — > c must be computed with non-zero charm mass, 
and this results in a complicated expression. [ ^0) In the 
SACOT scheme, the charm mass can be set to zero so 
that the final result for these sub-processes reduces to 
the very simple massless result. 

While the SACOT scheme allows us to simplify the 
calculation, the obvious question is: does this simpli- 
fied version contain the full dynamics of the process. 
To answer this quantitatively, we compare prediction 
for F$ harm obtained with 1) the SACOT scheme at or- 
der a\ with 2) the predictions obtained with the origi- 
nal ACOT scheme, 3) the ZM-VFN procedure in which 
the charm quark can appear as a parton but has zero 
mass, and 4) the FFN procedure in which the charm 
quark has its proper mass but does not appear as a 
parton. For simplicity, we take (j, = Q. 

In Fig. |^ we show F%{x, Q) as a function of Q for x = 
0.1 and x = 0.001 using the CTEQ4M parton distribu- 
tions. [ |H II We observe that the ACOT and SACOT 
schemes are effectively identical throughout the kine- 
matic range. There is a slight difference observed in 
the threshold region, but this is small in comparison 
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to the renormalization/factorization \x- variation (not 
shown). Hence the difference between the ACOT and 
SACOT results is of no physical consequence. The 
fact that the ACOT and SACOT match extremely well 
throughout the full kinematic range provides explicit 
numerical verification that the SACOT scheme fully 
contains the physics. 

Although we have used the example of heavy quark 
leptoproduction, let us comment briefly on the impli- 
cations of this scheme for the more complex case of 
hadroproduction.[ [lj ^3], |24|, ^5) At present, we have 
calculations for the all the 0(a 2 s ) hadroproduction sub- 
processes such as gg — ► QQ and gQ — ► gQ. At O(af) 
we have the result for the gg — > gQQ sub-processes, 
but not the general result for gQ — > ggQ with non- 
zero heavy quark mass. With the SACOT scheme, we 
can set the heavy quark mass to zero in the gQ — > ggQ 
sub-process and thus make use of the simple result al- 
ready in the literature.^] This is just one example of 
how the SACOT has the practical advantage of allow- 
ing us to extend our calculations to higher orders in 
the perturbation theory. We now turn to the case of 
heavy quark production for hadron colliders. 

5. Heavy Quark Hadroproduction 



There has been notable progress in the area of 
hadroproduction of heavy quarks. The original NLO 
calculations of the gg — > bb subprocess were performed 
by Nason, Dawson, and Ellis [ |23|], and by Beenakker, 
Kuijf, van Neerven, Meng, Schuler, and Smith[ p4| . 
Recently, Cacciari and Greco [ ^6) have used a NLO 
fragmentation formalism to resum the heavy quark 
contributions in the limit of large pt', the result is a 
decreased renormalization/factorization scale variation 
in the large pr region. The ACOT scheme was applied 
to the hadroproduction case by Olness, Scalise, and 
Tung.[ ^5) More recently, the NLO fragmentation for- 
malism of Cacciari and Greco has been merged with 
the massive FFN calculation of Nason, Dawson, and 
Ellis by Cacciari, Greco, and Nason, [ [27]]; the result 
is a calculation which matches the FFN calculation at 
low pr, and takes advantage of the NLO fragmentation 
formalism in the high px region, thus yielding good be- 
havior throughout the full pj- range. This is displayed 
in Fig. H where we see that this Fixed-Order Next- 
to-Leading-Log (FONLL) calculation displays reduced 
scale variation in the large pr region, and matches on 
the the massive NLO calculation in the small pr re- 
gion. Further details can be found in the report of the 
LHC Workshop b-production group^\ 

6. W + Heavy Quark Production 



> 



w 1000 — 



> 
o 



F0 

FONLL (MS) - 



3 

ft 

T3 



b 



500 - 



100 — 



50 




20 



40 60 

Pt (GeV) 



80 



100 



Figure 3. Differential cross section for b-production 
vs. pt comparing the Fixed-Order (FO) and the 
Fixed-Order Next-to-Leading-Log (FONLL) result in 
the MS scheme. The bands are obtained by vary- 
ing independently the renormalization and factoriza- 
tion scales. The cross section is scaled by m T with 
niT = y/ml +p T , and *fs = 1800 GeV, nib — 5 GeV, 
y = 0, with CTEQ3M PDF's. Figure taken from Cac- 



ciari, Greco, and Naso n, Ref. [ 27 



s For a related idea, 
Cacciari and Greco [ 



j the fragmentation function formalism of 
I in the following section. 



PDF Set 


Mass (GeV) 


LO 


WQQ 


NLO 


CTEQ1M 


m c =1.7 


96 


20 


161 


MRSD0' 


m c =1.7 


81 


20 


138 


CTEQ3M 


m c =1.7 


83 


20 


141 


CTEQ3M 


nib— 5-0 


0.17 


9.09 


9.33 



Table 1 

The W + charm-tagged one-jet inclusive cross section 
in pb for LO, W+QQ, and NLO (including the W+QQ 
contribution) using different sets of parton distribution 
functions. Table is taken from Ref. fEal. 



The precise measurement of W plus heavy quark 
(W+Q) events provides an important information on 
a variety of issues. Measurement of W+Q allows us 
to test NLO QCD theory at high scales and investi- 
gate questions about rcsummation and heavy quark 
PDF's. For example, if sufficient statistics are avail- 
able, W +charm final states can be used to extract 



9 The LHC Workshop b-production group is organized by 
Paolo Nason, Giovanni Ridolfi, Olivier Schneider, Giuseppe 
Tartarclli, Vikas Pratibha, and the report is currently in prepa- 
ration. The webnap-c for the h-nrodur.tion group is located at 
http://home.cern.en/n/nason/www/lhc99/ 
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Figure 4. Differential da/dplp for 7 plus tagged heavy 
quark production as compared with Pythia and the 
NLO QCD results. Figure taken from Ref. [ f29|. 

information about the strange quark distribution. In 
an analogous manner, the W+bottom final states are 
sensitive to the charm PDF; furthermore, W+bottom 
can fake Higgs events^ and are also an important back- 
ground for sbottom (b) searches. 

The cross sections for W plus tagged heavy quark jet 
were computed in Ref. [ , and are shown in Table. |l|. 
Note that this process has a large Zf-factor, and hence 
comparison between data and theory will provide dis- 
cerning test of the NLO QCD theory. While the small 
cross sections of these channels hindered analysis in 
Run I, the increased luminosity in Run II can make 
this a discriminating tool. For example, Run I pro- 
vided minimal statistics on W+Q, but there was data 
in the analogous neutral current channel 7+Q. Fig. || 
displays preliminary Tevatron data from Run I and the 
comparison with both the PYTHIA Monte Carlo and 
the NLO QCD calculations; again, note the large K- 
factor. If similar results are attainable in the charged 
current channel at Run II, this would be revealing. 

Extensive analysis the W+Q production channels 
were performed in Working Group I: "QCD tools for 
heavy flavors and new physics searches," and we can 
make use of these results to estimate the precision to 
which the strange quark distribution can be extracted. 
We display Fig. | (taken from the WGI report [ ^) 
which shows the distribution in x of the s-quarks which 
contribute to the W+c process.^ This figure indicates 

10 For a detailed analysis of this work including selection crite- 



Figure 5. Distribution of Events /0.01 vs. x of the 
s-quarks which contribute to the s + W — ► c process. 
Figure taken from Ref. [ |3(| . 

that there will good statistics in an x-range compara- 
ble to that investigated by neutrino DIS experiments; [ 
||, |j| hence, comparison with this data should provide 
an important test of the strange quark sea and the 
underlying mechanisms for computing such processes. 

7. The Strange Quark Distribution 

A primary uncertainty for W +charm production 
discussed above comes from the strange sea PDF, s(x), 
which has been the subject of controversy for sometime 
now. One possibility is that new analysis of present 
data will resolve this situation prior to Run II, and 
provide precise distributions as an input the the Teva- 
tron data analysis. The converse would be that this 
situation remains unresolved, in which case new data 
from Run II may help to finally solve this puzzle. 

The strange distribution is directly measured by 
dimuon production in neutrino-nucleon scattering .p| 
The basic sub-process is vN — > u~cX with a subse- 
quent charm decay c — ► [i + X' . 

The strange distribution can also be extracted indi- 
rectly using a combination of charged (VF 1 * 1 ) and neu- 
tral (7) current data; however, the systematic uncer- 
tainties involved in this procedure make an accurate 
determination difficult^ |3l| The basic idea is to use 

ria, see the report of Working Group I: "QCD Tools For Heavy 
Flavors And New Physics Searches," as well as Ref. [ p0[ . 
11 Presently, there are a number of LO analyses, and one NLO 
analysis. [ g g] 
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Figure 6. Variation of x s(x, n) for three choices of /j,, 
and also with a "SR" (slow-rescaling) type correction: 
x -> x(l + ml/Q 2 ). 
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the relation 

F 2 NC _ 5_ f 3 (s + s) - (c+ c) + 
F 2 CC ^ 



IS 



(1) 



to extract the strange distribution. This method is 
complicated by a number of issues including the xF% 
component which can play a crucial role in the small- 
x region — precisely the region where there has been a 
long-standing discrepancy. 

The structure functions are defined in terms of the 
neutrino-nucleon cross section via: 

d 2 a v,v G%me r , 9 «, 

— = ^-[^-^+^±^(1-1: 

It is instructive to recall the simple leading-order cor- 
respondence between the F's and the PDF's:p| 



F, 
xF. 



(v,D)N _ 



{u + 

= x {u — 



d + d + 2s + 2c} 
d-J±2s=F2c} 



(2) 



Therefore, the combination AxF^: 



AxF 3 = xF» N - xF* N = 4x{s - c} (3) 

12 To exhibit the basic structure, the above is taken the limit of 
4 quarks, a symmetric sea, and a vanishing Cabibbo angle. Of 
course, the actual analysis takes into account the full structure.] 



Figure 7. AxF 3 /2 vs. Q 2 for three choices of x. 
Calculations provided by S. Kretzer. 



can be used to probe the strange sea distribution, 
and to understand heavy quark (charm) production. 
This information, together with the exclusive dimuon 
events, may provide a more precise determination of 
the strange quark sea. 

To gauge the dependence of AxF% upon various fac- 
tors, we first consider xs(x, /i) in Fig. ^, and then the 
full NLO AxF 3 in Fig. [?]; this allows us to see the 
connection between AxF^ and xs(x, /i) beyond leading 
order. In Fig. we have plotted the quantity xs(x, fi) 
vs. Q 2 for two choices of £ in a range relevant to the 
the dimuon measurements. We use three choices of the 



fj, 2 scale: {Q 2 ,Q 2 



,P% }. The choices Q 2 and 



differ only at lower values of Q ; the choice 



Q 

but lies above for x 

12 i iiCD" „ c ,„o ,,2 



is comparable to Q and Q 



i? at x = 0.08 



0.015. The fourth curve labeled 
Q 2 + U SR" uses /j, 2 = Q 2 with a "slow-rescaling" type 
of correction which (crudely) includes mass effects by 
shifting x to x(l + m 2 /Q 2 ); note, the result of this 
correction is significant at large x and low Q 2 . 

In Fig. we have plotted the quantity AxFz/2 for 
an isoscalar target computed to order a\. We dis- 
play three calculations for three different x-bins rele- 
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vant to strange sea measurement. 1) A 3-flavor cal- 
culation using the GRV98[ |§ distributions ^ and 
H = \/ Q 2 + m 2 . 2) A 3-flavor calculation using the 
CTEQ4HQ distributions, and fJ, = Q. 3) A 4-flavor 
calculation using the CTEQ4HQ distributions, and 

H = Q- 

The two CTEQ curves show the effect of the charm 
distribution, and the GRV curve shows the effect of 
using a different PDF set. Recall that the GRV calcu- 
lation corresponds to a FFN scheme. 

The pair of curves using the CTEQ4HQ distribu- 
tions nicely illustrates how the charm distribution 
c(x,/x 2 ) evolves as \n(Q 2 /m 2 ) for increasing Q 2 ; note, 
c(x,/j, 2 ) enters with a negative sign so that the 4- 
flavor result is below the 3-flavor curve. The choice 
fi = Q ensures the 3- and 4-flavor calculation coin- 
cide at fi = Q = m c ; while this choice is useful for 
instructive purposes, a more practical choice might be 

fi ~ V<2 2 + ™ 2 : C /-> SeC - i and Ref - [O- 

For comparison, we also display preliminary data 
from the CCFR analysis. [|3l|| While there is much free- 
dom in the theoretical calculation, the difference be- 
tween these calculations and the data at low Q values 
warrants further investigation. 

8. Conclusions and Outlook 

A detailed understanding of heavy quark produc- 
tion and heavy quark PDF's at the Tevatron Run II 
will require analysis of fixed-target and HERA data 
as well as Run I results. Comprehensive analysis of 
the combined data set can provide incisive tests of the 
theoretical methods in an unexplored regime, and en- 
able precise predictions that will facilitate new particle 
searches in a variety of channels. This document serves 
as a progress report, and work on these topics will con- 
tinue in preparation for the Tevatron Run II. 

This work is supported by the U.S. Department 
of Energy, the National Science Foundation, and the 
Lightner-Sams Foundation. 
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